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ABSTRACT 

We present timing solutions and spin properties of the young pulsar PSR B0540—69 from analysis 
of 15.8 years of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent 
timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also 
perform fully coherent timing over large subsets of the data set in order to arrive at a more precise 
solution. In addition to the previously reported first glitch undergone by this pulsar, we find a 
second glitch, which occurred at MJD 52927 ± 4, with fractional changes in spin frequency Xv jv = 

(1.64 ±0.05) X 10“® and spin-down rate Xi/jv = (0.930 ±0.011) x 10“^ (taken from our fully coherent 
analysis). We measure a braking index that is consistent over the entire data span, with a mean value 
n = 2.129 ± 0.012, from our partially coherent timing analysis. We also investigated the emission 
behavior of this pulsar, and have found no evidence for significant flux changes, flares, burst-type 
activity, or pulse profile shape variations. While there is strong evidence for the much-touted similarity 
of PSR B0540—69 to the Crab pulsar, they nevertheless differ in several aspects, including glitch 
activity, where PSR B0540—69 can be said to resemble certain other very young pulsars. It seems 
clear that the specific processes governing the formation, evolution, and interiors of this population 
of recently born neutron stars can vary significantly, as reflected in their observed properties. 

Subject headings: pulsars: general — pulsars: individual (PSR B0540—69) — stars: evolution 


1. INTRODUCTION 

Several neutron stars (NSs) are presently observable 
as pulsars relatively soon after the supernova event as¬ 
sociated with their formation. These objects, which 
include well studied sources such as the Crab pulsar 
(PSR B0531±21) and the Vela pulsar (PSR B0833-45), 
appear to have been born with relatively high rotation 
frequencies, and undergo a rapid period of spin-down 
early in their lives as NSs. We can generalize the rate at 
which this rotational energy loss occurs via a power-law 
model: 

V = ( 1 ) 

where v and v are the pulsar spin frequency and fre¬ 
quency derivative, respectively, K is generally related to 
the NS moment of inertia, misalignment between the spin 
and magnetic axes, and magnetic field strength, and n 
is the braking index. For the generally utilized case of a 
pure magnetic dipole radiation mod el for pulsar spin- 
down, n = 3 for constant K (e.g., Ostriker & Gunn 


19691, but alternative models of the pulsar magneto¬ 
sphere, crust, and/or interior may predict different values 
for n, including a magnetic field that evolves with time 


(e.g., Blandford fc Romani|1988 Manchester & Peterson 


I989|), vary i ng spi n-magnetic inclination angle (e.g.,|Lyne 

eFar|[M3l[M5t , or other mechanism s cont ributing to 
angular momentum loss (Menou et al. 20011, such that 
K in Equation [2 is no loiiger assumed to be constant. 

Most pulsars nave such a small value of second fre¬ 
quency derivative that its determination is usually not 
feasible, even over very long timescales. In contrast, 
young pulsars are often observed to spin down at a suffi¬ 
ciently fast rate from their birth spin rates that the sec¬ 
ond derivative of the rotation frequency ii is measurable. 

^ rferdman@physics.mcgill.ca 


In these cases, we may take the derivative of Equation 
and find that we may cast the braking index as follows: 


n = 


( 2 ) 


We can therefore determine the braking index n through 
measurements of the rotation frequency and its deriva¬ 
tives. There are currently nine pulsars for which the 
braking inde x has been measured: all have been found to 
have n < 3 (Lyne et al.| 1993 1996 Middleditch et al. 


20061 [Livingstone et al. |2U(J7i |2(ji iT Espinoza et al.|2011a| 


Weltevrede et al.||20111 Roy et al.|20i2p . This indicates 

that there may indeed be a more complex picture regard¬ 
ing the magnetic field behavior and/or early evolution of 
the NS. Gontinued measurement of pulsar braking in¬ 
dices is therefore a key component to understanding NS 
magnetospheric behavior. 

There exist effects intrinsic to the NS than can influ¬ 
ence the reliability of the measurement of h, and thus n, 
consequently limiting the number of pulsars for which 
braking index can be measured. We observe timing 
noise as a long-term , quasi-random evo lution in pulsar 
spin frequency (e.g., Hobbs et al. 2010). Timing noise 
is thought to be related to magnetospheric instabilities, 
having been shown to have a s trong relationship to pro¬ 
file changes (Lyne et al. 2010). In contrast, glitches are 
near-instantaneous changes in rotation frequency, and of¬ 
ten in spin-down rate, which are believed to result from 
a sud den coupling of stellar crust and crustal superfluid 
(e.g., Espinoza et al.||2011b and references therein). As 
of this writing, there are 165 pulsars which have been ob¬ 
served to glitch, totaling over 470 recorded glitch events, 
according to the catalog maintain ed by Jodrell Bank Ob - 
servatory pulsar research grouip] ([Espinoza et al.|20lTb). 


http://www.jb.man.ac.uk/pulsar/glitches.html 
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It has been shown that the glitch event rate of pulsars 
is roughly correlated with age, such that younger pul¬ 
sars undergo more frequent glitch events t hat their older 
counterparts (e.g., Espinoza et al. 2011b). Indeed, very 
few glitches are seen from pulsars with large character¬ 
istic age; for example, only a single glitch has ever been 
repor ted from a millisecond pulsar (Cognard & Backer 
2004). The very youngest pulsars do, ho wever, appear to 


show a reduced glitch activity (see, e.g., Espinoza et al. 


2011b and references therein). One postulated cause tor 


this IS the relatively high temperatures in their interiors. 
In this scenario, after some cooling time, the pulsars will 


cent” 

Vela pulsar ( 

McKenna & Lyne 

1990 

Link et al. 

1999 

Dodson et al. 

5DI]2ir 


with a characteristic age = vj^lv ^ 1700 years, and a 
rotation frequency v = 19.8 Hz, or period Pspin = 50 ms. 
Its X- ray pulsations were discovered by [Seward et al.| 


(1984) in Einstein Observatory data associated with the 


Large Magellanic Cloud supernova remnant 0540—693. 
Its very large spin-down rate not only identified it as a 
likely young NS, but also made it possible to measure 
its braking ind ex. It was soon discovered to have opti¬ 


cal pulsations (Middleditch & Pennypacker 1985), and 
was immediately identified as very similar to the Crab 


Middleditch et al. 

1987 Hirayama et al.|2002|, the char- 

acteristics ot its o 

Dserved radio emission and giant pulse 

behavior (e.g., Manchester et al.||1993a 

Johnston &: Ro- 

mani|20U3 Johnston et al.|2UU4|), as wel 

1 as in the nebula 

Manchester et al. 

that its wind apparently powers (e.g., i 


~K braking index was first determined for 
PSR B0540—69 using both X-ray and optical mea¬ 
surements of the spi n frequency and its deri vatives, 
to be n = 3.6 ± 0.8 (Middleditch et al. 1987). This 
measurement disagreed sign ificantly with that of 


Manchester & Petersonj (|198^ , who subsequently found 
n ^ 2.01 ± 0.02 using independent observations of 
optical pulsations from PSR B0540—69. Since then, five 
additional braking index measurements have been made 
for this pulsar, few of which agree within measurement 
uncertainties (Middleditch et all |1987[ [Manchester 


fc Peterson 1989 Reiman fc Hasinger| 1990 Deeter 

et al.||19991 [Zhang et al.||2001[ jCusumaiio et al.||2003[ 

Livingstone et al. 2005a). The reason for this is likely 
due to a combination of timing noise contamination of 
the output parameter values, and possibly the existence 
of unaccounted-for glitch activity, particularly in the 
earlier X-ray and optic al data. Even in stud ies by 
Zhang et al. (|2001l) and jCusumano et al.j (|2003|), both 
of which performed phase-coherent timing of data from 
the Rossi X-ray Timing Explorer (RXTE) to derive 
pulsar parameters, neither resulted in consistent braking 
indices, nor did they agree on the cause of observed v 
and discontinuities near MJD 51325, with the former 
authors attributing it to a glitch, and the latter to 
strong timing noise. 

The most recent timing analysis of PSR B0540—69 by 
Livingstone e t al. ( 2005a ) r e-exam ined the RXTE data 
sets used by Izha ng et all (|2001|) (1.2 ye ars of RXTE 
observations) and Cusumano et al.|(2003) (1.2 years of 


RXTE data; ^ 8 years combined data set), and extended 
them to include 7.6 years of observations. In addition to 
a fully phase-coherent analysis, they performed partially 
phase-coherent timing, which divides the data set into 
shorter, non-overlapping subsets of data in order to de¬ 
termine the change in v and f' in a way that is more 
robust against the contamination of parameter determi¬ 
nation from timing noise. They confirmed tha t there 
is indeed a glitch at the same epoch claimed by jZhan^ 
et al. (2001). However, they report a significantly differ¬ 
ent nieas urement of the brak ing index of PSR B0540—69 
to that of IZhang et al.||2 001[ The fully coherent analysis 
performed byjLivingstone et al. (2005a I on the same data 
set supports the glitch parameter and loraking index mea¬ 
surements determined through their partially coherent 
analysis, and validating the reliability of that meth od. 

In this work , we re- analyze the data set used by [Liv-j 
ingstone et al.j (|2005a|) and extend its span to include all 
PSR B0540—69 data taken with RXTE until its decom¬ 
missioning, giving us a 15.8-year baseline for our timing 
work. In we describe the observations and our reduc¬ 
tion of the data set, including our time-of-arrival (TOA) 
calculations. In we present out timing efforts, and 
the results from both partially and fully coherent timing 
analyses, which we describe in more detail. We discuss 
our search for radiative changes , in cluding burst activ¬ 
ity in the RXTE time series in §4.1[ and attempt to set 
limits on emission variations from PSR B0540—69. Fi¬ 
nally, we devote to examining of the implications of 
our results, including a discussion of the much-published 
comparisons between PSR B0540—69 to the Crab pulsar, 
and how our findings might fit into the overall picture 
regarding the apparent similarity of these two systems. 
We also discuss how such comparisons may contribute to 
the understanding of the rotational and magnetic field 
evolution of PSR B0540—69, and possibly other NSs like 
it. 

2. OBSERVATIONS AND DATA REDUCTION 

All observational data of PSR B0540—69 used in this 
work were taken with the Proportional Counter Array 
(PCA; Jahoda et al. 20061) onboard RXTE. The PCA 
is an array of hve proportional counter units (PCUs), 
which together have a field of view of approximately 1 ° 
FWHM and a total collecting area of 6500 cm^. It is 
sensitive to an energy range of 2 — 60 keV, and records 
photon arrival times with a resolution of 1 /rs. We used 
those data acquired in “GoodXenon” mode, which si¬ 
multaneously provides high time resolution with energy 
information resolved over 256 spectral channels for each 
event that is not rejected due to background filtering. 

Correction of all event data to barycentric dynami¬ 
cal time (TDB) was applied using the barycorr script 
provided as part of the RXTE FTOOLS package^ 
Barycentering was perfor med relative to t h e puls ar co¬ 


ordinates determined by [Mignani et al.j (|2010|) with 
observations from the Hubble Space 1 'elescope, which 
remain the most precise position measurement for 
PSR B0540—69. For our analysis, we used events found 
in all PCU layers for our data set with energies within 

® Provided as part of the HEASOFT software analysis package; 
version 6.14 was used for this work (see http://heasarc.gsfc. 
nasa.gov/lheasoft/ftools/fhelp/barycorr.html) 
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the range 2 — 18 keV. 

Integrated pulse profiles were constructed for each ob¬ 
serving epoch by folding the given time series over 32 
phase bins. This was done by calculating the rotational 
phase of each event using the rotation frequency and fre- 
quency derivative reported by Livingstone et al. (2005a I 
as our reference ephemeris. We accumulated a prohle in 
this way for each observing epoch, and assigned a time 
stamp and pulse phase to the resulting profiles corre¬ 
sponding to the midpoint in time of the given epoch. 

To determine pulse times-of-arrival (TOAs) with which 
to perform our timing analysis, we first constructed a 
high signal-to-noise ratio (S/N) template prohle. This 
was done by aligning in rotational phase, then adding 
together all integrated prohles from our full data set. 
TOAs and corresponding uncertainties for a given epoch 
were then calculated by randomly generating 512 pro¬ 
hles based on the Poissonian count error on each phase 
bin for the corresponding integrated prohle. We then 
cross-correlated each trial prohle against the high S/N 
template in the Fourier domain, in order to calculate a 
phase shift relative to the fundamental harmonic of the 
template prohle. The mean shift from all trials for a 
given epoch was then converted to absolute time offset 
via the reference pulse period, and added to the time- 
stamp of eac h prohle to produce a corresponding TOA 
([Taylor 19921. The corresponding uncertainty for a given 
I'UA was found from the standard deviation of all cor¬ 
responding trial TOAs used for that observing epoch. 
Prohle smearing effects on later-epoch TOA uncertainties 
from folding event data with our reference ephemeris are 
negligible; the fractional rotation frequency error 8v jv is 
orders of magnitude smaller than the precision required 
to avoid signihcant spin phase drift within a prohle bin, 
over the typical integration times of < 5000 s. Uncertain¬ 
ties therefore rehect the diherences in prohle S/N, which 
is principally due to the different integration times, as 
well as the number of active PCUs, throughout the data 
set. 

3. TIMING ANALYSIS 

In th is work, we have extended the [Livingstone et ah] 


(2005a I data set by more than a factor of two, to include 
5.8 years of RXTE data. We have performed a phase- 
coherent timing analysis on these data, which allows us 
to account for every rotation of the NS. PSR B0540—69 
is an isolated pulsar; phase coherence can be achieved 
by htting the measured pulse TOAs to a spin model that 
determines the rotational phase (j) of the pulsar at a given 
time t, via a Taylor expansion of the spin frequency of 
the pulsar and its derivatives: 

cj){t) = (j){tQ)+VQ{t — tQ) + -I>Q{t — to)^ + -l'o{t—to)^ + . ■ ■ 

( 3 ) 

Here, to is a reference epoch in our model, and t'o, f'o, and 
i>o represent the pulsar spin frequency and its derivatives 
at epoch to- 

For all timing analyses, we use the tempo2 s oftware 
package ([Edwards et al.[|2006| [Hobbs et al.[[2006[) , which 
fits our bary centered 'i'UAs to a model describing the 
spin evolution of the pulsar. As with our pulse profile 
construction, we use d the spin frequency and fre quency 
derivatives from the Livingstone et al. (2005a I model 


ephemeris as a starting point for our timing analyses. 
We have mod i fied t he pulsar position to match that of 
Mignani et al. (2010| mentioned above, and have updated 
the Solar System model to use the JPL DE421 ephemeris 
( Standish|2004[) in order to account for the motion of the 
Earth. We have not included any parameters describing 
a priori knowledge of glitch behavior of this pulsar, al¬ 
lowing for an independent analysis and comparisons with 
previous results. 

Following the general procedure used by [Livingstone[ 
et al. (2005a), we performed timing analysis in two inde¬ 
pendent ways: in §3.1[ we describe our effort to reduce 
the effects of intrinsic pulsar spin noise by performing a 
partially coherent analysis, in which we divi de t he data 
set into relatively short timing intervals; in §3.2[ we de¬ 
scribe our fully coherent timing analysis, in which we 
attempt to phase-connect the entire data set (over two 
large overlapping sections), and compare our results from 
both methods. 

3.1. Partially Phase-coherent Analysis 

Young pulsars often exhibit high levels of timing noise 
that can inhibit accurate measurement of the intrinsic 
frequency derivatives, and thus braking indices, as well 
as identification and characterization of glitches. While 
phase-coherent timing can provide very precise parame¬ 
ter determinations, higher-order spin derivatives are of¬ 
ten needed to remove the effects of spin noise. This is 
typically a red process, occurring over long timescales, 
leading to covariances between these parameters. 

In order to mitigate the influence of timing noise, we 
divided the full data set into 59 non-overlapping sub¬ 
sets. The length of each subset was determined based on 
its timing properties; specifically, we included the num¬ 
ber of TOAs that allowed for a phase-coherent timing 
fit that required only the rotational frequency v and its 
first derivative v to maintain white, Gaussian-distributed 
post-fit timing residuals. For each subset, we have taken 
the center of the time span as the reference epoch for the 
timing fit. 

We find evidence for a glitch at MJD 5134 8 ± 36, con- 
sistent with that fo und by Livingstone et al. (2005a) and 
Zhang et al. (2001). In this analysis, we quote the mid¬ 
point between the epochs just before and after the dis¬ 
continuity found in i' as the value of the glitch epoch. 
The (symmetric) uncertainty is simply calculated to be 
the distance in time between the glitch and either neigh¬ 
boring epochs. In our extended data set, we have also 
found a second glitch at MJD 52925 ±54. In Figurej^we 
plot our measurements of the spin frequency derivative 
at each subset epoch for this partially coherent analysis; 
both glitches are denoted by a vertical dotted line at their 
respective epochs, and are the focus of the inset plots. Fi¬ 
nally, we find some evidence for a previously unreported 
glitch near the start of the data set, which cannot be cor¬ 
roborated due to a lack of earlier data points. If this is 
a true glitch, we determine that it would have occurred 
at MJD 50264 ± 68. 

We were able to fit three distinct slopes to our mea¬ 
surements of !>, allowing us to determine v for three sep¬ 
arate date ranges: before the first glitch (not including 
the first epoch); between our two confirmed glitch events; 
and after the second glitch. We determined uncertain¬ 
ties for our i) measurements through a bootstrap tech- 
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Fig. 1.— Partially coherent timing measurements of the spin frequency derivative v of PSR B0540—69. Points show fit values of u for 
each data subset. Vertical dashed green lines show the glitch epochs at MJDs 51349 it 36 and 52925 it 54. Blue dashed lines show fit slopes 
to !> data before the first glitch, between glitches, and after the second glitch. Insets focus on a portion of the data near each glitch, with 
red shaded areas showing the extent of the possible glitch epoch. 


nique. Here, we randomized the data order many times 
to construct a distribution of fit values, the width of 
which was taken to be the uncertainty in i). This analysis 
showed the errors taken directly from the output covari¬ 
ance matrices from our linear fits to be underestimated, 
and therefore report our bootstrap-derived uncertainties 
in Table [TJ Our measurements of v are consistent within 
their mutual Icr uncertainties; we therefore claim no sig¬ 
nificant third frequency derivative measurement can be 
made from our analysis of this data set. 

With our measured !>, and ii, we were then able 
to calculate a mean braking index for each epoch range, 
finding n = 2.127 ± 0.017, 2.127 ± 0.013, 2.135 ± 0.004 
for the first, second and third data spans, respectively. 
These are consistent over the data set within Icr; we 
therefore quote a mean value n = 2 .129 ± 0.012, also 
consiste nt with the reported value by [Livingstone et al.| 
(2005a I. For each glitch, we calculate the change in fre- 
quency derivative Aic by determining the difference in 
the extrapolated v from the partially coherent solutions 
before and after the glitch epoch. From these, we find 
a fractional change jv = (1.554 ± 0.006) x 10“^ at 
the first glit ch epoch, which is c onsiste nt with the value 
reported by Livingstone et al. (2005a). At the second 
glitch epoch, we find Ai>/i> = (9.460 ± 0.011) x 10“®. In 
Figure we plot a summary of our partially coherent 
timing analysis. The top two panels show the residual 
values for v and after their respective trends prior to 


the first glitch have been subtracted out. In the case of ic, 
we note that one can easily see the measurement discon¬ 
tinuities near the reported glitch epochs (including the 
proposed, but unconfirmed, initial glitch). We cannot, 
however, make a significant measurement of the intrinsic 
changes in rotation frequency Ai/ due to their small mag¬ 
nitudes. This is because the overall measured frequency 
change for both glitches is dominated by the change in 
frequency derivative, over and above the frequency step 
itself, on a timescale that is significantly shorter than the 
time between the epochs surrounding each glitch. The 
bottom two panels of Figure!^ show the measured values 
of V and corresponding braking indices for each of the 
three epochs, which are delimited by the two reported 
glitches. 

3.2. Fully Phase-coherent Analysis 

We performed a coherent timing analysis of our entire 
data set, in order to confirm the results of our partially 
coherent analysis. More importantly, it allowed us to ob¬ 
tain more precise measurements of the glitch epoch, and 
to determine the change in the pulsar rotation frequency. 
Ah' (and therefore derive the fractional change Av/v) at 
each glitch epoch, which we were not able to achieve in 
our partially coherent analysis. 

As stated in the previous subsection, long-term tim¬ 
ing of young pulsars is often affected by a great deal 
of red spin noise, and this is certainly the case for 
PSR B0540—69. In fact, it was not possible to phase- 
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TABLE 1 

Timing parameters for PSR B0540-69 


Data set 

Data span (yr). 

Date range (MJD). 

Number of TOAs. 

Solar system ephemeris model . 


15.8 

50123.2 - 55898.7 
2247 

DE421 


Set Quantities 

H. A.^ (J2000), a . 

Decl.‘‘ (J2000), 5. 


05:40:11.202 

-69:19:54.17 


Measured Quantities—Partially Coherent Timing i^’it 


Subset 1 

Subset 2 

Subset 3 


Date range (MJD) 50296.9 - 51333.1 51342.9 - 52921.0 52930.6 - 55830.9 

Second derivative of rotation frequency, i/ (xl0“^^ s“®) 3.79(3) 3.78(2) 3.773(6) 


Braking index, n. 2.127(18) 2.127(13) 2.135(4) 


Glitch epoch (MJD). 51349(36) 52925(54) 

Ai> (xlO-i"^ s-2). -2.9(2) -1.77(15) 

(xlQ-^)..1.554(6)0.9460(11) 


Measured Quantities—Coherent Timing I'it 



Subset 1 

Subset 2 

Date range (MJD) 

50123.2 - 52935.7 

51342.9 - 55898.7 

Rotation frequency, v (s“^). 

19.8024438758(9) 

19.7746860321(7) 

First derivative of rotation frequency, v (xl0~^^ . 

-1.8780249(6) 

-1.8727175(8) 

Second derivative of rotation frequency, v (xl0~^^ s~^) 

3.721(10) 

3.772(5) 

Reference timing epoch (MJD). 

51197 

52910 

Braking index, n . 

2.08(3) 

2.12(3) 

RMS timing residual (ms). 

0.74 

1.75 

of fit. 

1293.22 

10088.45 

Degrees of freedom. 

719 

1432 

Number of fit frequency derivatives. 

12 

12 

Glitch epoch (MJD). 

51335(13) 

52927(4) 

Ai/ (xlO-* s-i). 

2.5(6) 

3.24(10) 

Avjv (x 10“®). 

1.3(3) 

1.64(5) 

Ai> (xl0“^^ s“®). 

-2.75(7) 

-1.74(2) 

Ai>/i> (xlO-"^). 

1.46(4) 

0.930(11) 


^ Adopted from [Mignani et al.| ||2010[| 


connect the entire data set with an RMS of the timing 
residuals within a pulse period. This was due to the 
need for more than 12 frequency derivatives to do so, 
which would exceed machine representation capabilities 
of those values, and is thus not available in Tempo 2. We 
thus divided the data set into two overlapping segments: 
the first subset includes all data up to, but not including, 
the second glitch, and the second subset includes all data 
after, but not including, the first reported glitch. 

Figure shows the post-fit timing residuals for both 
subsets ofthe full data set, with an artificial vertical off¬ 
set added between the two sets for clarity. In both cases 
we have included several high-order frequency deriva¬ 
tives in our ht in an attempt to “whiten out” as much 
as possible the effects of red timing noise; specifically, 
we employed 12 derivatives for both subsets. A sum¬ 
mary of our coherent analysis results can be found in 
Table We report second frequency derivatives from 
each subset that are inconsistent by an amount greater 
than their mutual 3(t uncertainties; however, our use of 
high-order frequency derivatives due to timing noise ef¬ 
fects has likely introduced covariances between spin pa¬ 
rameters that would affect our i) measurements. This 
also casts doubt on coherent-based determination of a 
third derivative v that clearly represents its long-term 
spin evolution, rather than red-noise effects. 

Our analysis resulted in a braking index measurement 
for each subset, evaluated at each corresponding refer¬ 
ence epoch, also reported in Table Simply adopting 


the measurement errors based on the formal fit uncer¬ 
tainties would likely underestimate the resulting prop¬ 
agated uncertainty on the braking index n; this is due 
to timing noise contamination that results in covariances 
between higher-order frequency derivatives. We there¬ 


fore follow the more conservative method taken by Liv¬ 


ingstone et al. (2005a), by taking the uncertainties on 
each n to be the standard deviation of the calculated n 
resulting from a set of timing fits that use between 2 and 
12 frequency derivatives. This better reflects the sensi¬ 
tivity of n to the number of parameters in our fit. We find 
n = 2.08 ± 0.03 and 2.12 ± 0.03 for the first and second 
subsets, respectively, in agreement with each other, and 
with the results of our partially coherent timing within 
Icr uncertainties. This consistency between subsets and 
analyses gives us confidence in the unchanging nature 
of the braking index over the nearly 16 years span of 
this data set. We therefore quote a mean braking index 
n = 2.10 ± 0.03 from this fully coherent timing analysis. 

We were additionally able to confirm both glitch 
epochs found with our partially coherent analysis, and 
measure the fractional changes in both rotation fre¬ 
quency and spin-down rate through coherent timing anal¬ 
ysis. As expected, these glitches have resulted in low- 
level changes in i/; we find jv = 1.3 ± 0.3 x 10“® and 
1.64 ± 0.05 X 10“® for the first and second glitches, re¬ 
spectively; we also find Ai>/i> = 1.46 ± 0.04 x 10“"^ and 
0.930 ± 0.011 X 10“"^. While our determined fractional 
change in the rotation frequency agrees well within Icr 
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Epoch(MJD) 


Fig. 2. — Summary of partially coherent timing analysis of 
PSR B0540—69. Panel (a): residual measurements of rotation fre¬ 
quency, determined by removing our reference model, which is valid 
prior to the first reported glitch (fiducial epoch MJD 51197). Dis¬ 
continuities are seen to occur around the time of each glitch. Panel 
(b): residual spin-down rate where the fit slope prior to the first 
glitch is first removed. Clear offsets are seen to be associated with 
the reported glitches, as well as persistent spin-down rate increases 
after each glitch. Panels (c) and (d): second rotation frequency 
derivative and braking index, respectively. The values (dashed 
line) and uncertainties (dotted lines) shown were determined us¬ 
ing data taken before the first glitch, between the two reported 
glitches, and after the second glitch; these date ranges represented 
as colored regions. 


with the reported value of Livingstone et al. (2005a I, this 
is not quite the case for the measured fractional spin- 
down increase, only agreeing near to their mutual 2cr 
ranges; however, our spin-down glitch sizes Ah' do agree 
with those determined from our own partially coherent 
solutions. This is not entirely surprising, given the high 
level of timing noise seen in the post-fit residuals from 
the coherent timing analysis (Figure^, which results in 
glitch parameter values that vary with the number of 
spin frequency derivatives included in our timing model. 
We did not find a glitch at the beginning of this data set, 
near MJD 50264; this may be due to both a lack of data 
around that epoch and the effects of timing noise, or else 
because the glitch did not in fact occur. 


4. SEARCHING FOR RADIATIVE CHANGES 

We conducted an analysis of our data set to deter¬ 
mine whether there existed any flux variations, and if 
so, whether these could be associated with the observed 
glitches. To do so, we first filtered our event data in order 


to exclude known high-background epochs that may have 
otherwise biased our flux measurements. This included 
the rejection of data during South Atlantic Anomaly pas¬ 
sage, Earth occultation and bright Earth effects, and/or 
electron contamination. We also rejected epochs with 
pointing offsets greater than 0.02°. Finally, we did not 
use any data taken with PCUO on or after 2000 May 
12 (MJD 51676), or with PCUl on or after 2006 De¬ 
cember 25 (MJD 54094) due to loss of their respective 


propane lay ers ^see, e.g., Janoaa et ai. zIUUd vorcoran 
|et al.|[2gT0] as well as the KATE guest observer online 
facilitjQ; while this does not affect our TOA calculations 
and timing analysis, it does influence the calibrated flux 
values, causing the interpretation of PCUO and PCUl 
flux data to be unreliable after the above dates. 


4.1. Flux Measurements 

We calibrated our flux data separately for each PCU 
based on the collimator response measured for the Crab 
pulsar as a function of angular separation between the 
pulsar positi on and the nominal telescope pointing (as 
measured by |Jahoda et al.||2006|) at each given epoch. 
In calculating huxes, we determined the pulsed count 
rate for eac h PCU using the RMS flux variations. This 
method (see An et al(]|2013 Archibald et al.|[2014 for a 
more detailed description) results in flux measurements 
with lower bias than the usual method of integrating the 
pulse profile above the background levels; the latter re¬ 
lies on the uncertainty in the perceived minimum of the 
profile, which results in larger biases for noisier profiles. 
In Figure]^ we show the measured RMS pulsed flux over 
the data set for PCU2. We find that there are no signif¬ 
icant flux enhancements that are clearly associated with 
the glitch events reported in 


4.2. Profile shape changes 

In order to test for changes in profile morphology, we 
used the template profile we created for our timing anal¬ 
ysis described in m as a model for the profile shape of 
PSR B0540—69. We first accumulated pulse profiles from 
data taken within each of the subse ts us ed for the par¬ 
tially coherent analysis described in §3.1| This was done 
in order to obtain high S/N profiles, while allowing us to 
have a reasonably sized sample of resulting pulse profiles 
over the time span of this data set. For each subset pro¬ 
file, we scaled it to match the template after removing 
the mean value of both profiles, in order to account for 
any relative flux changes. We then aligned each profile 
in pulse phase with the template, through a frequency- 
domain cross-correlation similar to that outlined in 0 
for TOA calculation. 

We determined the goodness-of-fit of our template pro¬ 
file to each aligned and scaled added profile by calculat¬ 
ing the statistic. In each case, we find values that are 
consistent with a distribution with the relevant num¬ 
ber of degrees of freedom; this includes profiles represent¬ 
ing epochs just before and after glitches, which have fits 
that yield p-values greater than 0.99. We therefore detect 
no significant change in profile shape over our entire data 
site, including at those epochs related to the observed 
glitches. Based on the lack of measurable flux or pro¬ 
file shape changes, and using the measured glitch epochs 


^ http://heasarc.gsfc.nasa.gov/docs/xte/ 
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Fig. 3.— Fully phase-coherent timing residuals for PSR B0540—69, divided into two overlapping subsets. The first of these includes all 
data before the second reported glitch, and the second subset comprises all data after the first glitch. Scales for the first and second timing 
analysis subsets are shown on the left and right vertical axis labels, respectively. Glitch epochs are denoted by vertical dashed lines. 


determined from coherent timing, we find the maximum 
time for radiative changes, and subsequent decay, result¬ 
ing from a glitch to be 21 and 13 days for the first and 
second glitches, respectively. 

4.3. Search for Burst Activity 

We also conducted a search for bursts over the nearly 
16 years for which we have PSR B0540—69 data. In order 
to do so, we created a time series from the event data us¬ 
ing 1/32-s time bins, in the same 2 — ISkeV energy range 
with which our timing analysis was conducted. Assum¬ 
ing a Poisson distribution for our time-series data, we 
looked for individual data points that were outliers from 
running means over 10, 100, and 1000 s subsets of our 
time series, using in cremental steps of 50% of each of 
these timescales (see Gavriil et al.||2004 Scholz & Kaspi 
|201H where descriptions of our employed method are 
given). We did not find any significant burst activity in 
this data set. 

The lack of X-ray flux variation, burst activity, or pro¬ 
file changes in PSR B0540—69 is consistent with the over¬ 
all observed behavior of most rotation-powered pulsars, 
which have stable fluxes and pulse shapes. Exceptions to 


this are PSR J1846—0258 (Gavriil et al. 2008 Kuiper & 
Hermsen||2009[ |Livingstone et al.||2010| Archibald et al. 
2015t andTSR J1119-6 127 dWeltevrede et al.||20TT 


Antonopoulou et al. ]M5t , which have an approximately 


order-of-magnitude larger dipolar magnetic helds than 
PSR B0540—69), and have displayed radiative changes 
near a glitch. This suggests a possible link to magnetars, 
which also ofte n show flux variations connected to ob - 
served glitches (Dib et al.|2007 j2008| Dib & Kaspi|2014 ). 
The young, hig 


i-held pulsar PSR 1 


509—58 has shown 


no gl itch activity in over 28 years (Livingstone & Kaspi 


20111, however, if one is revealed in future observations 


it could help to shed important light on the properties of 
young pulsars with high magnetic fields relative to other 
observed young NSs. 

5. DISCUSSION 
5.1. Braking Index 


Our partially and fully coherent analyses of 
PSR B0540—69 result in consistent determination 
of braking indices, which appear to remain constant 
throughout our 15.8-yr data set, within measured 
uncertainties (see Table [^. Furthermore, our measured 
values of braking index are consistent within Icr of 
the quoted values of n = 2.149 ± 0.009 reported by 
Livingstone et a l. (2005a), and n = 2.125 ± 0.001 by 
Gusumano et al. 1 20031; our values are not consistent at 
the m utual 3(7 level with that reported by |Zhang et al.| 
(20011, which was likely contaminated by the effects of 


timing noise. 

As with all published braking index values, that of 
PSR B0540—69 is significantly lower than the value of 
3 predicted for an ideal magnetic dipole held, indicating 
a departure from the standard model of dipole radiation 
for pulsar slow-down. Unlike studies of the Crab pul¬ 


sar and PSR B1509—58 , for example (Livingstone et al. 
2005b|[L^et al.|2015D , we have unfortunately not been 


able to measure a reliab le th ird frequ ency derivative v 
for PSR B0540—69 (see |3.1|and |3.2|), which could have 
given further insight into the possible time variation of its 
magnetic held. However, recent modeling of young pul¬ 
sar braking indices shows that the known braking indices 
may be explained by higher-order held evolution due to 
Hall drift in the NS crust, causing an enhancement of 


the d ipole held component (Gourgouliatos & Gumming 
2015|). Further discoveries of pulsars with measurable 


braking indices are therefore crucial for dehning a better 
picture of the magnetic held evolution of young pulsars. 

5.2. Glitches or timing noise? 

There are several factors that give us conhdence that 
our reported events are in deed glitches: (a) our conh rma- 
tion of the glitch found by Livingstone et al. (2005a I near 


MJD 51348, as well as the consistency in its correspond¬ 
ing measured fractional change in v, lends credence to 
its existence, rather than attr ibuting it to timing nois e 
effects, as was postulated by Gusumano et al. (2003|); 
(b) our discovery of a second glitch near iVlJD 52925 of 
similar fractional change in both i/ and to that of the 













































































































Ferdman et al. 


first glitch demonstrates consistent behavior; and (c) the 
persistent offsets observed in spin-down rate after each 
glitch leads us to favor the interpretation of these discon¬ 
tinuities as genuine glitches. 

We can compare this behavior to that of the Crab 
and other young pulsars. Although some have experi¬ 
enced higher glitch event rates, some have undergone 
many low-level glitches, several of which have smaller 
fractional change in both v and v than those seen in 
PSR B0540—69. Figure plots lS.v/v and Af'/f' for 
PSR B0540— 69, as well as those compiled by Espinoz^ 
et al. (2011b| for the Crab and other pulsars. The glitch 
sizes we hna in the PSR B0540—69 data do not seem 
to occupy a particularly unique part of parameter space, 
and are well within those glitch event sizes found for 
many young pulsars. Finally, as discussed in the fol¬ 
lowing section, the Crab pulsar is also known to show 
persistent steps in frequency derivative. 


5.3. A Crab twin? 

Since the discovery of PSR B0540—69, much has been 


discussed of its resemblance to the Crab pulsar (Seward 
et al. 1984[), and their similarity in several properties 


has supported this comparison. For the Crab, its char¬ 
acteristic age Tc,Crab = 1258 years, measured spin lumi¬ 
nosity Ecrab = 4.5 X 10^®ergs“^, and surface magnetic 
field strength Rg.Crab = 3.8 x 10^^ G, are indeed tan- 
talizingly close to those of PSR B0540—69, which has 
A ,0540 = 1672years, Eo 54 o = 1-5 x 10^®ergs“^, and 
Ss,o 540 = 5-0 X 10^^ G, hinting at a similar formation and 
evolution. Additionally, as discussed above, the Crab has 
experienced many low-level glitches events in both n and 
z>, several of which are even smaller in fractional size than 
those observed in PSR B0540—69. The Crab has, how¬ 
ever, undergone several large glitch events, which have 
not yet been seen in PSR B0540—69. While this dif¬ 
ference may be attributable to a variety of NS interior 
processes, it should be noted that the Crab pulsar has 
been observed with nearly four times as long a time base¬ 
line, and with much denser and more regular cadence, 
than has PSR B0540—69. It is also therefore possible 
that PSR B0540—69 has experienced unobserved large 
glitches, and/or may do so in the future. It should be 
noted that a 36% increase has recently been observed in 
the spin down rate of PSR B0540—69. This is thought 
to be due to a state change perhaps caused by magneto- 
spheric processes, rather than a glitch, since no accom¬ 


panyi ng rotation-rate change was seen (Marshall et al. 


2015 

While the Crab has many more observed glitches 
than PSR B0540—69, with 24 observed ov er 45 years 


2011b Lyne et al. 2015), both pul- 
low the relatively low obs erved glitch 


(Espinoza et al. 
sars appear to to. 
rate of th e very youngest pulsar s (e.g., Shemar & 
Lyne 1996 Espinoza et al. 2011b). Along with the 
Crab and PSR B0540—69, this group also includes 


PSR J1119-6127 {Tr ^ 1700 years: Camilo et al. 2000 


Weltevrede et al.|2011 Yu et al.|20l3 

), which has shown 
ions, and the sim- 
(tc ~ 1700 years; 

only 3 glitches in 16 years of observal 

ilarly young pulsar, PSR B1509—58 

Kaspi et al.||1994| |Livingstone et al.|f 

2005b[ |Livingstone| 

& Kaspi||20iip, which has not been o 

Dserved to undergo 


Year 

1996 1998 2000 2002 2004 2006 2008 2010 



MJD 


Fig. 4.— RMS pulsed fluxes for observations of PSR B0540—69, 
for the energy range 2 — 18keV; shown are flux data from PCU2. 
Green dashed lines show glitch epochs, and orange dotted lines 
show the epochs at which PCUO (MJD 51676) and PCUl (MJD 
54094) lost their respective propane layers. 


low glitch activity in PSR B0540—69 may therefore not 
be all that surprisin g, and may be related to its high in¬ 


ternal temperature (McKenna & Lyne 


1990), compared 


to NSs with roughly an order of magnitude higher char¬ 
acteristic age, such as the Vela pulsar, which has had 
signihcantly larger glitches. However, the relatively high 
glitch activity of magnetars, which appear to be the very 
hottest of all NS s, is not consistent with this hypothesis 
(Dib et al. 2008). Finally, it may also be possible that 
PSR B0540—69 underwent other low-level glitch events 
that, due to high-amplitude timing noise , were impossi- 


ble to detect in earlier timing efforts (see Espinoza et al. 
|2014[ for a thorough discussion of small glitches and tim- 
ing noise). However, as seen in Figure^ the glitches we 
observe for PSR B0540—69 have fractional sizes that are 
consistent with the glitch size distribution of the Crab 
pulsar. 

We can perform a rough comparison of glitch size and 
frequency between the Crab pulsar and PSR B0540—69 
through the so-called glitch activity parameter Ag = 
y)(Ai//i/)/At, which quantifies the cumulative fractional 



possibility of an unobserved period of frequent and/or 
heavy glitch activity for PSR B0540—69 outside the dates 
observed in this work—that the Crab pulsar glitches are 
more frequent, and larger in magnitude, than those of 
PSR B0540—69, and that this is one property in which 
these two objects differ signihcantly. 

As can be seen in Figure ^ ib measurements in our 
partially coherent analysis show persistent changes af¬ 
ter both glitches. This is also a feature of Crab pulsar 


glitches (Gullahorn et al.||1977 

Demiahski & Proszyhski 

1983 

Lyne et al. 

20T^ 

), which also show persistent in- 


plete recovery that terminates before reaching pre-glitch 
z> values. Recoveries are not seen in our PSR B0540—69 
data, probably due to a much lower density data set, 
compared to the impressi ve observing cade nce obtained 
for the Crab over 45 years ( Lyne et al.|2015 ). The persis¬ 
tent steps in i> have been explained by models in which 
the effective NS moment of inertia is decreased due to 
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Fig. 5.— Fractional frequency and frequency derivative glitch 
sizes found for PSR B0540—69 in this work (blue stars) com- 
pared with those found for the Crab pulsar (black filled circles 
|Lyne et al.|[20T^ , as wel l as all glitches (gray') reported in the 
glitch catalog published by|Espinoza et al.|(|2011b| (and maintained 
online at http://www.jb.mairracTuE7puTsar7gTitches.html) that 
fall within the above plot limits. 


the trapping of vortices in the interior, which are only 
released when a glitch occurs (Alpar et al. 1994). An 
alternative mode l invokes the po st-glitch rnovement of 
crustal “plates” (|Ruderman|| 1991 1) , to which the surface 
magnetic field is coupled, causes an increased misalign¬ 
ment of the magnetic dipole moment from the spin axis, 
thus increasing the e xternal torque a nd producing the 
observed offsets in v ( Link et al.||1992 ). 

Finally, as with all measured braking indices, both the 
Crab pulsar and PSR B0540—69 have n that is signifi¬ 
cantly less than 3, indicating a departure from the stan¬ 
dard picture of magnetic dipole braking. Along with a 
contribution from magnetospheric wind torques, one pos¬ 
sible explanation involves the secular change in pulse pro¬ 
file component separation observed in the Crab pulsar, 
indicating a movement of the magnetic dipole axis away 
from th e pulsar rotation axis, aff ecting its spin-down evo¬ 
lution (Lyne et al. |2013 2015). However, the paucity 
of detected X-ray photon events for PSR B0540—69 re¬ 
quires us to choose a much coarser pulse profile binning 
in order to maintain sufficient S/N for timing analysis 
and emission studies. This makes it nearly impossible to 
discern the subtle changes in pulse morphology that have 
been observed in the radio Crab pulsar prohles. Further¬ 
more, the weak r adio flux of PSR B0540—69 ([Manch¬ 
ester et al.|1993a ) makes it difficult to perform this type 


of analysis at those wavelengths, as has been done with 
the Crab pulsar. However, future radio telescopes and 
instrumentation with high projec ted sensitivity, su ch as 
the Square Kilometer Array (e.g., [Watts et al.|2015|) may 
provide data sets that allow for such a study. 


6. CONCLUSIONS 

We have extended the timing analysis of the very young 
pulsar PSR B0540—69 to include a time span of nearly 16 
years of RXTE data. We have confirmed the epoch and 
properties of the first g litch, previously reported by [Liv-[ 
ingstone et al. (2005a), and discovered a second glitch. 
We have conhrmed that the braking index of this pul¬ 
sar is consistent in value throughout our data set, and 
that no significant radiative changes have observed to be 
associated with its glitch activity, a behavior that is in 
line with the majority of young pulsars. Many of the ob¬ 
served properties of PSR B0540—69, such as fractional 
glitch depth, are consistent with several of those seen in 
the Crab pulsar. However, we do find that quantities de¬ 
scribing other aspects of its behavior, such as the glitch 
rate, size distribution, and activity parameter, can dif¬ 
fer substantially from the Crab, while showing greater 
similarity to those of other young pulsars. In any case, 
there may be a variety of processes that contribute to the 
spin behavior of PSR B0540—69 and other recently born 
NSs. Further monitoring of these sources and discov¬ 
ery of more young pulsars will be crucial toward gaining 
a better understanding of their evolution and observed 
properties in the context of the overall NS population. 
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